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Abstract 

A recent study of gamma-ray data from the Galactic Center motivates the investigation of light 
(~ 7 — 10 GeV) particle dark matter models featuring tau lepton pairs as dominant annihilation 
final state. The Lepton-Specific two-Higgs-doublet Model (2HDM-L) provides a natural frame- 
work where light, singlet scalar dark matter can pair-annihilate dominantly into tau leptons. We 
calculate the nucleon-dark matter cross section for singlet scalar dark matter within the 2HDM-L 
framework, and compare with recent results from direct detection experiments. We study how 
direct dark matter searches can be used to constrain the dark matter interpretation of gamma ray 
observations, for different dominant annihilation final states. We show that models exist with the 
correct thermal relic abundance that could fit the claimed gamma-ray excess from the Galactic 
Center region and have direct detection cross sections of the order of what needed to interpret 
recent anomalous events reported by direct detection experiments. 
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I. INTRODUCTION, MOTIVATION AND OUTLINE 



These are exciting times in the search for the particle nature of dark matter. A remark- 
able convergence of experimental capabilities is connecting direct detection campaigns with 
indirect searches and with the collider program (see e.g. Ref. [I]). Tantalizing signatures, 
and a few controversial claims, have focused some attention on a specific class of weakly 
interacting massive particle (WIMP) dark matter candidates, with relatively light masses 
(roughly, in the range between 7 and 10 GeV). Here, we briefly review these experimental 
data, and we give a detailed proof of existence of a particle setup (originally proposed in 
Ref. j2]) that can successfully explain them, while at the same time having a thermal relic 
density matching the observed dark matter density in the universe. 

After reporting an irreducible excess of low-energy (below 3 keVee) events in a p-type 
point contact Germanium detector [3J, the CoGeNT Collaboration recently released the 
results of fifteen months of cumulative data-taking. These data provide evidence for the 
presence of a modulated component, with a statistical significance of 2.8a, of unknown ori- 
gin, compatible in principle with a light-mass WIMP [3], featuring a mass around or slightly 
below 10 GeV, and a spin-independent scattering cross section with nucleons on the order 
of 10 -40 cm -2 . The WIMP interpretation is in tension with results from the CDMS-II ex- 
periment, that recently carried out a dedicated low recoil-energy threshold analysis [5], and 
with results from the XENON10 [7j and XENON100 \9\ experiments. Ref. [10], how- 
ever, critically argues that the limits presented in [7] and in [9] are, for various and different 
reasons, overstated, and that the parameters space compatible with a WIMP interpretation 
of the CoGeNT signal is presently not ruled out by either result. Notice that this claim 
is at present highly controversial - Ref. [IS] and [IB] , for example, evaluate and essentially 
refute the possibility of a low L e ^ envisioned in Ref. JTU] as a caveat to the XENON results 
ruling out the CoGeNT signal. In addition, the note [H] argues that due to entirely differ- 
ent background cuts for the two detectors, the CDMS-II and CoGeNT measured irreducible 
spectra might actually be in agreement, and compatible with a ~ 7 GeV WIMP signal. 
Remarkably, the CoGeNT signal is also in principle compatible with the annual modulation 
observed by DAMA/LIBRA if interpreted as WIMP recoils [12] (see also [13]). Clearly, 
further experimental information is in order to settle the question, including closer inspec- 
tion of detector issues, and new data from other experiments, such as CRESST, that also, 
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tentatively, presented tantalizing signal excesses [14] . 

Interestingly, an excess gamma-ray emission from the general direction of the Galactic 
Center region has been claimed recently to be present in the first two years of data collected 
by the Fermi-Large Area Telescope (LAT) |15j . Under the assumptions of Ref. [15] for the 
diffuse and point-source background, the morphology and spectrum of the emission is in 
principle compatible with a 7-10 GeV mass WIMP annihilating dominantly into a pair of 
tau leptons, with a thermally averaged cross section times relative velocity in the range 
(av) ~ 5 x 10 -27 to 5 x 10 -26 cm 3 /s. Ref. [16] also recently compared the Fermi LAT data 
with predictions within the constrained minimal supersymmetric extension of the Standard 
Model. 

Ref. [17] also recently entertained the intriguing possibility that the origin of the claimed 
gamma-ray excess from the center of the Galaxy might be related (via the produced rela- 
tivists electrons and positrons synchrotron-radiating at radio frequencies in the magnetic 
fields of the Galaxy) to an explanation for the intensity, spectrum and morphology of the 
so-called WMAP haze - the unexplained excess at radio frequencies correlating with the 
innermost few degrees from the Galactic Center [18] . 

It must be noted that data from the Galactic Center region appear, from the analyses 
and arguments presented in Ref. [T$l |2"U] to be consistent with a purely astrophysical origin, 
possibly connected with physics in the region surrounding the central supermassive black 
hole Sag A*. In particular, Fermi LAT gamma-ray data from the Galactic Center region 
was found in [21] to be entirely compatible with a standard background model consisting 
of a diffuse component plus point sources with different spectra from what was assumed 
in Ref. [15] . As is the case for the direct detection results, a more in-depth and detailed 
analysis and modeling of the gamma-ray data in the complex and crowded region of the 
Galactic Center seems timely and in order. 

Given the tantalizing claims for WIMP signals, it is by all means tempting to accommo- 
date those signals within a coherent particle-physics framework for particle dark matter. A 
model-independent analysis of particle physics scenarios that could accommodate the results 
reviewed above was recently given in Ref. [22] . A specific, minimal realization of a beyond- 
the Standard Model dark matter framework that also could closely match the claimed direct 
and indirect detection excesses was outlined in Ref. [2] (for a closely related supersymmetric 
model see also Ref. [23]). 
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Here, we are concerned with a closer and detailed inspection of the "lepton-specific" two- 
Higgs-doublet model with singlet scalar dark matter proposed in Ref. [2]. Following [2J, we 
will hereafter indicate this scenario with the acronym 2HDM-L. This class of models was, 
in a broader context, previously considered in Ref. [21], with an emphasis on larger masses 
and on a dark matter interpretation of the charged leptonic cosmic ray anomalies reported 
by PAMELA and ATIC Early work on the phenomenological implications of general two- 
Higgs doublet models, including a classification of the relevant Yukawa interaction structures, 
was given in Ref. [3E]- More recently, the phenomenology of 2HDM with general Yukawa 
structures have been addressed in several studies, including Ref. [39H4T] . A recent study 
that addressed the dark matter phenomenology of a class of models closely related to the 
one under investigation here was given in Ref. [12] . The 2HDM-L plus singlet structure can 
also lead to phenomenologically viable baryogenesis and neutrino mass generation models, 
as discussed in Ref. [15] . 

As we review in the following section, the model consists of a scalar singlet plus two 
SU(2) doublet scalar fields, one coupling to the quark sector and one coupling to the lepton 
sector. Ref. [2j uses approximate expressions for the pair-annihilation cross section and for 
the relic density, and argues that matching the claimed gamma-ray excess and the observed 
dark matter relic density via thermal production might be possible. Again resorting to ap- 
proximate expressions, Ref. [2] argues that the direct detection cross section for the 2HDM-L 
dark matter model with a dark matter particle around 7-10 GeV predicts a spin-independent 
dark matter-nucleon cross section which is one order of magnitude below the required cross 
section to explain the CoGeNT anomalous signal. 

Ref. [2] studies in some detail the collider signatures of the 2HDM-L model, showing that 
the Standard-Model-like Higgs generically acquires a large invisible width to the dark matter 
scalar, while the lepton-friendly Higgs decays predominantly invisibly with a few percent 
branching ratio into t + t~. While discovery of the predicted (heavier) Higgs structure is 
challenging at the LHC, Ref. [2] argues that it would be straightforward at an e + e~ linear 
collider with a 500 GeV center of mass energy. 

In the present study, we carry out a comprehensive study of the 2HDM-L parameter 
space, including detailed and complete expressions for the relevant cross sections and for 
the dark matter thermal relic abundance. Our main result is that we find several models 
with the right mass, spin- independent and pair-annihilation cross sections to explain the 
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direct and indirect signals reviewed above, and with a thermal relic density matching the 
observed dark matter abundance in the universe. We therefore argue that the 2HDM-L 
setup is particularly compelling as a framework that would give a comprehensive WIMP 
interpretation to recent experimental and observational claims. 

The remainder of this study consists of a review of the 2HDM-L model including con- 
straints from electroweak precision measurements and perturbativity (sec. Inj), details on the 



relevant direct and indirect cross sections (sec. Ill), our model parameter space scan results 
(sec. 



IV), and, finally, our discussion and conclusions (sec. M). 



II. THE MODEL 

The model considered in [2] contains two complex SU(2)-doublet scalar fields Hl and 
Hq that, through a Yukawa mechanism, generate the masses of the charged leptons and 
the quarks, respectively, and a gauge-singlet real scalar field S, whose corresponding mass 
eigen-state is the dark matter candidate. A Z 2 symmetry is imposed on the SU(2) fields, 
under which only Hl and em (right-handed leptons) are odd. A distinct unbroken global Z 2 
symmetry under which S is odd is also assumed to ensure the stability of the dark matter 
candidate. The scalar potential of the model is given by : 



V = HxH^Hq + fi 2 HlH L + ^S 2 + XqS 2 HqHq + \ l S 2 H{H l 
+ \ S S" + \ x {H{Hl) 2 + Hrf Q H Q ) 2 + X 3 (H^H Q )(HlH L ) 
+ ~^~((HqHl) 2 + h.c.) + X 5 (HqHl H l Hq). (1) 
After electroweak symmetry breaking and the minimization of the potential we can write: 

//, I ^ .... A>*Q=( . H ' 2+ ..^^\- (2) 



+ H[ + iA'J/y/2 / ' \(v 2 + H 2 + iA' 2 )/V2 



Assuming CP conservation, the couplings Ai_5 and the vacuum expectation values t>i_2 are 
taken to be real. In what follows, un-primed particles denote mass eigenstates. 

A. The Mass spectrum 

The minimization of the tree-level potential leads to the following mass spectrum: 
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M 2 Hi = v 2 Xx + v 2 2 X 2 - ^Hvlvl + ivlXi-vlXzY; (3) 

M 2 H2 = v\X x + v 2 2 X 2 + sJl 2 v\v\ + [v\X x -v\X 2 ) 2 \ (4) 

M\ = -v 2 X 4 ; (5) 

M\ t = -^ 2 (A 4 + A 5 ); (6) 

M 2 = ^(2m 3 + vlX L + v 2 2 X Q ), (7) 

where v 2 = v \ + v\ = (246 GeV) 2 and L = A 3 + A 4 + A 5 . The mixing in the neutral CP-even 
sector is defined as: 

Mh x = cos a H[ + sin a H 2 , 

M H2 = — sin a H[ + cos a H' 2 . (8) 

We parametrize the model with the following set of independent variables: 

M Hl , M H2 , M A2 , M H ±, M s , X L , A Q , tmP = v 2 /v u a. (9) 

The additional singlet quartic parameter is set to A s = 0.1 throughout our analysis, since 
the ensuing phenomenology is widely independent of it. Additionally, for the purpose of 
the phenomenological study we are pursuing, the masses of the pseudo-scalar and charged 
states play an insignificant role and will be fixed to 120 GeV for the benchmark study. We 
are thus left with 6 free parameters. We refer the Reader to Ref. [2] for additional details 
on the model, including Feynman rules for the couplings of the scalar sector of the theory 
to quarks, leptons and gauge bosons. 

B. Constraints on the model 

In addition to the requirement that the scalar potential be bounded from below, we 
imposed the following theoretical and experimental constraints on the model: 

• Electroweak precision tests 

The oblique parameters S, T, U provide tight constraints on generic theories beyond 
the Standard Model [25]. For models with additional scalars only, the T parameter 



6 



is the most relevant one, as scalars offer negligible contributions to S and Uf2El |2~T] . 
The stringent bounds on the T parameter from electroweak measurements |28j are: 

-0.08 <T< 0.14. 

Following the notation of [29], the T oblique parameter for the 2HDML-L model is : 

AP = 64/2^2 [f(M hP M A2 ) + 3(F(M W , M Hl ) - F(M Z , M Hl )) 

cos(a - /3) 2 (F(M^ + , M H2 ) + F(M A2 , M H2 ) + 3F(M W} M Hl ) - 3F(M Z , M H% )) 
sin(a - (3) 2 (F(M H +, M Hl ) + F(M A2 ,M Hl ) + 3F(M W , M H2 ) - 3F(M Z , M Ha )) 
= a em T, (10) 

where a em is the fine-structure constant and the function F is defined as (x,y > 0): 

{x + y xy x 
— In - x ± y, 
2 x-y y (n) 
<^> x = y. 

Perturbativity 

Following Ref. [10] we impose perturbativity bounds on all the couplings of the model: 

A l <7<4vr i = l,...,5 (12) 

This requirement leads to an upper bound on the masses of the higgses of ~ 850 
GeV. The lower bound is conservatively taken to be 114 GeV, in accord with LEP-II 
searches for the Standard Model Higgs [30]. The upper bound on the masses of the 
charged and pseudo-scalar states is ~ 630 GeV. For a discussion of parameter space 
bounds from perturbative unitarity, see also Ref. 



Flavor Physics 

Extremely significant constraints on this model come from flavor physics, and in par- 
ticular from measurements of BR(6 — > 57), which in the presence of a charged Higgs 
boson receives important new physics contributions. We used here the results of 
Ref. [40J to implement that constraint. We find that other flavor observables are less 
constraining for this model than BR(6 — > 37). 
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III. DARK MATTER IN THE 2HDM-L MODEL 



As detailed in the introduction, recently Hooper and Goodenough [15] argued that an 
excess gamma-ray emission might exist in the inner region of the Galactic Center (within 
approximately 1.25°) and be due to a dark matter particle in the mass range 7—10 GeV 
annihilating primarily to tau leptons. We start our analysis by focusing on the slice with 
Mg = 8 GeV of the parameter space. We then randomly scan over the other free parameters 
in the model (Eq. (|9| above) and apply the constraints on thermal relic density, indirect de- 
tection cross-section and the latest direct detection spin-independent cross-section exclusion 
bounds. We make no approximations in the expressions for all relevant cross sections. 

As regards dark matter detection channel other than direct detection under scrutiny here, 
it is worthwhile to mention neutrino telescopes, that in principle provide very stringent 
constraints to models such as those under investigation here. Super-Kamiokande data can 
provide rather tight bounds on models with the mass range and pair- annihilation final states 
we considered here, as recently shown in Ref. [37]. While some of the models considered 
here might be in tension with Super-Kamiokande data, it is difficult to discern it directly 
from e.g. Fig. 3 in Ref. [37], given that for the present models one has an admixture of 
s- and p-wave annihilation, and a mixture of final states that have significantly different 
constraints from neutrino telescopes. 

For a particle to be a viable cosmological dark matter candidate - barring episodes of 
entropy injection (see e.g. [31J) - the pair annihilation cross section must be large enough 
to prevent over-production in the early Universe via standard thermal freeze-out. For the 
particular mass range we are interested in, and considering the Yukawa suppression for light 
fermions, the main kinematically open channels are annihilations to tau leptons and to b 
quarks via Hi and H 2 s-channels. Notice that there are no co-annihilation channels open in 
this model. 

In order to study the viable region of this benchmark we perform a random scan over 
the 6-dimensional parameter space, fixing, as specified above, Ma 2 = M H + = 120 GeV. We 
compute the dark matter relic abundance and the spin-independent scattering cross section 
off of protons with the micrOMEGAs code [32, 33j - Notice that to maximize the number of 
points passing the constraint on the T parameter we however tune, when necessary, M H + 



to get a value in the correct range (in practice, we solve Eq. (10) finding the value of M H -t 
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FIG. 1: Spin-independent DM scattering cross section off protons versus annihilation cross-section 
times the relative velocity of the DM particles. Vertical lines guide the eye towards models with 
good cosmological relic density, whereas horizontal ones represent the "CoGeNT region" and the 
lower bound of XENON100. Regions where r + r~ final state accounts for less than 50% and 
more than 60% and 80% of the total annihilation channels are shown. We fix Mg = 8GeV and 
M A2 = M H + = 120 GeV. 

for which Ap = 0, for values of M H + compatible with the LEPII constraints on the charged 
Higgs boson mass). Again, the tuning of this parameter does not affect quantities relevant 
for dark matter phenomenology. 

In Fig. [I] we show the result of a scan in the plane [asr{crv)). Depending on the pa- 
rameters of the model, the final states of dark matter annihilation (via s-channel Hi and 
H 2 exchange) are mainly tau leptons and b-quark pairs. Notice that in the absence of reso- 
nances or coannihilation, a thermal relic density corresponding to the observed dark matter 
density is found for models with a pair-annihilation cross section of a few xlO -26 cm 3 /s. To 
guide the eye, we indicate the relevant range with vertical lines in the plot. 
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The fit to gamma-ray data obtained in Ref. [TH] requires at least 60% of the final states 
to be tau leptons. We explicitly show with different colors points where t + t~ final state 
accounts for less than 50% (red) and more than 60% (blue) and 80% (green) of the total 
annihilation channels. Currently, the experiment that, at face value, puts the most stringent 
exclusion limits is XENON100 [9] that reaches a spin-independent cross-section of ~ 4 x 
10~ 43 cm 2 for a dark matter particle of 8GeV. It is clear from the plot that the 2HDM-L 
model has large parameter space regions compatible with cosmology (relic density) and with 
the XENON100 null results, while featuring a dominant t + t~ annihilation final state. 

Interestingly enough, we find regions that are compatible with the annual modulation 
signal reported by the CoGeNT collaboration [I], although imposing the relic density con- 
straint reduces significantly the region with a r + r~ annihilation final state dominating with 
more than ~ 60% branching ratio. The fact that the b-quark-dominated region is a thin 
line (i.e. there exists a strong correlation between direct and indirect dark matter detection 
rates), whereas the tau-dominated regions are wider on the plane is due to the different 
contributions entering the two precesses of annihilation and the scattering cross-section off 
protons. Indeed, for the scattering off protons the parameter dependence of the vertices 
(except the mass of the quarks) is essentially the same as for the annihilation into quarks, 
since it is basically the same process in a t-channel. We have then that the ratio a si/ (erf) 
is constant. This argument does not hold when the annihilation proceeds through leptons 
because, in that case, the two SU(2) Higgses couple differently to leptons and quarks. 

An illustrative analytical derivation of this result is as follows: The annihilation cross- 
sections into fermions can be approximated as [2] 



{av) 



N c m 



2 r /1™2 n 3/2 



Aix 



1 i - 



Cl (13) 



s 

where N c is the number of colors of the species f, sis the square of the center-of-mass energy 
and Cf includes all the contributions of the 6 parameters of the model and are given by: 
( Xq tan (3 cos a — Xl sin e*) sin a (Xq tan (3 sin a + Xl cos e*) cos a 



a 



(Xq cos a — Xl cot f3 sin cv) cos a (Xq sin a + Xl cot j3 cos o;) sin a 



Cq s-M% 2 s-M* Hl ■ (14) 

The spin-independent cross-section of scattering off the proton is instead: 

4 

TfL 
10 



Model 


1,1 X2 l 




Ar 


An 


sin ct 


tan B 


(Jqt fcm^l 


fi/i 2 


(av) [10~ 26 cm 3 /s] ] 


A 


114.8 


177.1 


6.75 


-0.55 


0.638 


5.67 


6 79xl0~ 45 


0.12 


2.22 


B 


114.7 


270. 


6.87 


-0.21 


0.01 


5.79 


2.17xl0~ 43 


0.12 


2.24 


C 


117. 


163. 


-3.9 


1.01 


0.169 


7.1 


8.69xl0~ 41 


0.12 


2.14 


D 


114.6 


162.3 


-0.48 


0.87 


0.25 


6.62 


8.04xl0" 41 


0.091 


2.96 



TABLE I: Examples of set of parameters for the Ms = 8 GeV benchmark 



where m p is the proton mass, C q is the same function that appears in Eq. (13) evaluated 
at s = 0, and is a constant function of nuclear form factors, for more detail see e.g 
Ref. [51]. Therefore, one has an approximate linear relation between a si and (crv), the 
slope of which depends only on the mass of the DM particle. We will use this result in the 
next section where we perform a full scan including varying Ms- Table |T] presents sets of 



parameters that pass the constraints listed in |II B[ have the cosmological relic density and 
annihilate predominately to tau leptons. Models A and B are compatible with the most 
stringent bounds from XENON100, whereas C and D are compatible with CoGeNT for 
different values of tan (3 and low-mass Higgs. 



IV. FULL PARAMETER SPACE SCAN 



We proceed now to a more general scan, to probe the low mass region of this model more 
systematically. We follow the same strategy as for the study of the benchmark but we relax 
the constraints on the masses Ms and Ma 2 = M H + = 120 GeV. We vary linearly Ms from 
5 GeV to 15 GeV, and we sample randomly 10 7 points in the following ranges: 

• Scalar Masses: from 114 GeV up to the relevant pertubative limit (~ 600 GeV). 

• Charged and Pseudo-scalar particles: from 100 GeV up to their perturbative limit 
(~ 450 GeV ). 

• tan 0: from to 50. 

• sin a: from to 1. 
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FIG. 2: Spin-independent dark matter scattering cross section off protons versus annihilation cross- 
section times the relative velocity. Vertical lines bound models with good cosmological relic density 
whereas horizontal ones represent the CoGent region. Regions where t + t~ final state accounts for 
less than 50% and more than 60% and 80% of the total annihilation channels are shown. 



We apply the constraints described in sec. |II B to every parameter space point. The result 



of the scan are shown on the (ogi, (crv)) plane in Fig. [2] and on the (Ms, <Jsi) plane in Fig. [3j 
In Fig [3] we focus on direct detection only, and show the experimental sensitivities of 
DAMA, CoGeNT, XENON100 and CDMS @]. Notice that different analyses obtain different 
regions compatible, at some given confidence level, with the DAMA/LIBRA modulation 
signal. For example, the shape of the parameter space compatible with DAMA/LIBRA 
obtained in e.g. Ref. jU] and [35J differ in the allowed mass range, and this effect depends 
upon different assumption on the choice of the bins used, the quenching factors, the velocity 
distribution, and the definition of the confidence levels. Here, we choose to employ the 
results of Ref. [35], which emphasize the possibility that the CoGeNT and DAMA signals 
might have a common origin, which is an especially appealing scenario for the model we 
investigate here. Regions below the green line contain 2HDM-L models that are consistent 
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5 6 7 8 9 10 11 12 

M s [GeV] 



FIG. 3: Spin-independent DM scattering cross-section off protons as a function of the mass of DM 
with the regions of DAMA/LIBRA (without channeling) [35] and COGENT [I] and the sensitivity 
lines of CDMS and XENON100. For XENON100: (A) corresponds to the collaboration limits 
and (B) to a sensitivity obtained under different assumptions [10]. For CDMS: (A) is the CDMS 
limits [36] and (B) corresponds to CDMS-II results [5], disputed in Ref. [TO]. The blue lines delimit 
the region where lepton final states start dominating and the green one delimits the region where 
the tau lepton final state dominates with more than 80%. 

with having the correct thermal relic abundance, and featuring a pair-annihilation final state 
branching ratio larger than 80% into the t + t~ final state. Models below the dark blue line, 
instead, have branching ratios larger than 50%. The figure illustrates that we find ample 
parameter space in the 2HDM-L framework to accommodate the signals reported in direct 
dark matter detection with dark matter models with the right relic abundance that feature 
large branching ratios into t + t~. 

V. DISCUSSION AND CONCLUSIONS 

We investigated the phenomenology of the dark matter candidate arising from a lepton- 
specific two-Higgs doublet model with a gauge-singlet scalar. We imposed on the parameter 
space both electroweak constraints and limits from LEP searches for the Higgs. We then 
calculated the spin-independent dark matter scattering cross section off of nucleons, and 
the pair-annihilation cross section. The two quantities have a very simple correlation for 
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dominant b-quark annihilation final state, while the correlation is lost for tau-lepton final 
state. Using the micrOMEGAs code we calculated the thermal dark matter relic abundance. 

We reached here somewhat different, and more optimistic conclusions than what reported 
on a recent study on the same model [2]. Specifically, we showed here that within the 
model under consideration one can account for the anomalous signals reported by CoGeNT, 
DAM A/LIBRA and the claims of a gamma-ray excess from the Galactic Center, while, at 
the same time, having a thermal dark matter relic abundance compatible with the observed 
universal dark matter density. Should the various anomalous signals discussed here pass 
forthcoming scrutiny, the simple dark matter framework we investigated in this study might 
indeed deserve additional attention. 
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